Metal-insulator transitions in systems with electron-phonon and Coulomb interactions 
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We consider a model which includes the electron-phonon coupling to A s or Jahn- Teller H 9 
phonons, the Coulomb interaction U and an exchange integral K. We study the metal-insulator 
transition for the integer fillings n = 3 and 4. We find that the coupling to A fl phonons increases the 
critical value U c where the transition takes place, while the coupling to H 9 phonons decreases U c . 
Without electron-phonon coupling the Hund's rule coupling also decreases U c . There is, however, 
an interesting competition between the Hund's rule coupling and the Jahn- Teller effect. Thus the 
reduction of U c due to the H 9 phonons becomes smaller when the Hund's rule coupling is turned 
on. The implications of this is discussed in the context of A n C6o- 



I. INTRODUCTION 

The alkali-doped fullerenes, A„Ceo (A= K, Rb), raise 
many interesting problems due to the important role 
played by both the electron-electron and the electron- 
phonon interaction Q]. For n = 3 and n — 4 each al- 
kali atom is assumed to donate about one electron into 
a six-fold degenerate t± u band. A3C60 are metals and 
superconductors while A4C60 are insulators |jL||. 

The t\ u orbital is only partly filled for both systems, and 
band structure calculations predict both systems to be 
metallic ||. A4C60 must then be an insulator due to in- 
teraction effects neglected in band structure calculations. 
Under pressure A4C60 becomes metallic M, while some 
fullerenes at normal pressure and with the doping n = 3 
(NH 3 K 3 C6o, CS3C60) are not superconductors and prob- 
ably insulators, but become superconductors under pres- 
sure [p|-pO[| . This suggest that these systems are relatively 
close to a metal-insulator transition. A metal-insulator 
transition is usually discussed in terms of the ratio be- 
tween the Coulomb interaction U between two electrons 
on the same molecule and the one-particle band width 
W. For A3C60 and A4C60 this ratio is, however, almost 
identical [raljj]. It is then interesting to ask which fac- 
tors determine whether a system is on the metallic or 
insulating side. 

These systems have been studied in a Hubbard-like 
model, and it was found that the three-fold degeneracy 
of the tiu orbital plays an important role by increasing 
the ratio U/W where the metal-insulator transition takes 
place jl2|,[l3]]. It was furthermore found that the lattice 
structure is important for the difference between A3C60 
(fee) and A4C60 (bet) jljjl^]. Hubbard-like models pre- 
dict, however, that A4C60 is an anti-ferromagnetic insu- 
lator, while it is known experimentally that there are no 
moments in A4C60 [qjlq]. The electrons in A n Cgo have a 
relatively strong interaction with Jahn- Teller intramolec- 
ular phonons with H ff symmetry and a weaker interaction 
with A g phonons. The interaction with the Jahn- Teller 
phonons favor a low spin state and might lead to a non- 
magnetic insulator fT^ . At the same time there is, how- 
ever, a Hund's rule coupling, which favors a high spin 



state. This leads to an interesting competition between 
the Jahn- Teller effect and the Hund's rule coupling. The 
purpose of this paper is therefore to study the influence 
of the Jahn- Teller effect and the Hund's rule coupling on 
the metal-insulator transition. We do this in the context 
of the Fullerenes, but similar effects also occur in many 
other systems, e.g., transition metal compounds Jig] . 

In this paper we study a model of A„C6o which in- 
cludes the electron-phonon coupling and the electron- 
electron interaction. In A„Cgo partly occupied orbital t\ u 
couples to two nondegenerate phonons of A g symmetry 
and to eight five-fold degenerate (Jahn- Teller) phonons 
of H g symmetry. For A„Cgo the important coupling is 
believed to be to the H ff phonons and we therefore 
study a model with a coupling to one H 3 phonon. To see 
the effect of having a Jahn- Teller phonon, we also com- 
pare with a model containing an A g phonon. The models 
are solved in the dynamical mean-field theory (DMFT) 
|]l9f or by using exact diagonalization. To interprete the 
results we also study analytically a single molecule and a 
simple two-site model. A brief summary of some aspects 
of this work are published elsewhere p3 . 
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Ifl we present results for a free molecule and 



we discuss the typical parameter range. We 
present results for a two-site model in Sec. [fy| and de- 
scribe the DMFT calculation in Sec. fv[ T he results are 
given in Sec. [v| and discussed in Sec. |VI]| . 



II. ISOLATED MOLECULE 

We first study an isolated molecule. We consider the 
case of the coupling to an A g phonon and two cases 
of coupling to Jahn- Teller phonons, namely the E x E 
case, where a two-fold degenerate level couples to a two- 
fold degenerate phonon, and the T x H case, where a 
three-fold degenerate level couples to a five-fold degener- 
ate phonon. The E x E case is the simplest model with a 
Jahn- Teller effect, and the T x H case is a simple model 
of a Ceo molecule. Thus we consider the Hamiltonian 



H, 



cl — pt 



,t 



UJQ 



1 



+9 E E E E Kil' c ^ c -v(^ + (i) 



^ cr m 



where the first term describes the electronic level, the 
second the phonon and the third the electron-phonon 
coupling. The coupling matrices are determined by 
symmetry and they are given in, e.g., Ref. |2p|-|22]]. For 
the coupling to an A g phonon is diagonal with all 

the diagonal elements equal to unity, = 6 mm > . For 

the E x E case 



1 

-1 

and for the T x H case 
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The overall coupling strength is determined by g. In 
a metallic molecular solid with narrow bands, there is a 
simple relation between g and the dimensionless coupling 

a m 



Ah = |iV(0)^- 
6 loq 



for the T x H case and 



\ E = 2N{0)^- 



(4) 
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for the E x E case. Here iV(0) is the electron density of 
states per spin. We also define a Jahn- Teller energy 



Ejt — 1 — • 

UJ 



(6) 



In addition we include the Coulomb interaction 

Hjj = U xx 2J n m ^n m i + U xy na- m n a ' m ' 

m aa m<m 

+ E ^m^L'a'^rnv'^m'a ( 7 ) 

<j<j m^m 

where n ma = ^^ lo .'0 m(T is an occupation number, U xx and 
U xy are the Coulomb interactions between equal and un- 
equal orbitals, respectively and K is an exchange integral. 
The Coulomb integrals are related via 



U X y Uxx 2-K". 



(8) 



For A g phonons and a free molecule, the electron- 
phonon problem is reduced to a displaced oscillator, and 
we can easily obtain the electron-phonon contribution to 
the ground-state energy 



E Ag (N) = -N 2 E JT . 



(9) 



where E(N) is the energy of a free molecule with N elec- 
trons. The Jahn- Teller case is more complicated and can- 
not be solved exactly. We therefore focus on the lowest 
order contribution to E(N). For the Coulomb interac- 
tion this contribution is of first order in K and for the 
electron-phonon interaction it is of second order in g. 
The results are shown in Table | for the E x E case and 
m Table || for the T x H case. We have subtracted an 
average Coulomb energy, 

E(N) = E(N) - ^N(N - l)U av (k), (10) 

where U av (k) is the average interaction for the ground- 
state in a full shell with 2k electrons and the orbital de- 
generacy k. By definition, E(N) is then zero for a full 
shell. Results for the low spin state (singlet or doublet) 
and the high spin state (triplet or quartet) are shown in 
Table | 

Neglecting the Coulomb interaction (U and K) the re- 
sults for the T x H case can be written as |23j22| 



where 



E^(N) = -C(N)E JT , 



1 for TV = 1,5 
C(N) = \ 4 for N = 2,4 
3 for N = 3 



(11) 



(12) 



This is a weak-coupling result for the case when A (or g) 
is small. For the strong-coupling case |23|,|22|| 



Efj c (N) = C(N)E. 



JT- 



(13) 



In the strong-coupling case the prefactor in front of 
C(N)Ejt has been reduced by a factor of 5/2. 

From tables | and || it follows that the low spin state 
is lower in energy for K < 2Ejt in the E x E case and 
for K < \E JT in the T x H case. The results illustrate 
the competition between the Jahn- Teller effect and the 
Hund's rule coupling. According to Hund's first rule, it 
is favorable to form a high spin state to minimize the 
Coulomb energy. On the other hand, the Jahn- Teller 
effect favors a low spin state. For instance, if the first 
phonon (v = 1) of E symmetry (the left matrix in Eq. 
(^)) is excited, it is favorable to put two electrons (with 
opposite spins) in the level m = 2, resulting in a singlet. 
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The tables show that in the low spin state the energy 
E(N) is lowered by the electron-phonon coupling but is 
increased by the exchange integral K. In the high spin 
state, on the other hand, the exchange coupling lowers 
the energy, while the electron-phonon interaction lowers 
the energy less than in the low spin state or not at all. 

This can be further illustrated by considering the case 
of two electrons in the E x E system. For K < 2Ejt, 
the important states are 



\i) = 4Mi\ 



vac) 



|2) = v4 t v4j 



vac) 



(14) 



These states couple directly via the exchange interaction 
and indirectly via the electron-phonon interaction. Thus 
we also include states where one phonon has been excited. 
In the corresponding Hamiltonian matrix the part corre- 
sponding to the states with one phonon are "folded" into 
the part corresponding to i = 1 or 2, and we obtain 
the matrix 



U xx -6Ej T -2Ej T + K 
-2Ej T + K U XX -6E JT 



(15) 



If K < 2Ejt the corresponding energy is 



U xx - 6E JT + (-2Ej T + K) = U av {2) - 8E JT + -K, 

(16) 

which clearly shows the competition between the Jahn- 
Teller coupling and the Hund's rule coupling. For K > 
2Ejt a triplet state of the type t/'j+'^T l vac ) with the 
energy U av {2) — 4K/3 becomes the lowest state. 



III. PARAMETER RANGE 

It is now important to establish the parameter range 
appropriate for A„C6o- In particular, the relative size 
of g and K is important, since this determines whether 
the Jahn- Teller effect or the Hund's rule coupling wins. 
The exchange integral K has been estimated [|24fl from 
an ab initio SCF calculation |B5[. From the calculated 
multiplet splitting the result K — 0.11 eV was obtained. 
Correlation effects are expected to reduce this number, 
since correlation effects in particular lower the energies 
of the low spin states. The difference in Coulomb energy 
between the low spin states and the high spin states is 
then reduced and the effective K becomes correspond- 
ingly smaller. For instance, for some atomic multiplets, 
correlation effects were found to reduce the multiplet 
splitting by about 25 % j2(|. The electron-phonon cou- 
pling constants have been estimated from photoemission 
experiments for a free molecule (271] . Here we replace the 
eight H g modes by one effective mode. The frequency 
of this mode is chosen as the logarithmically averaged 
frequency 



A = \ v ; A lncjo = J]] A^ \r\u v 



(17) 



where A„ and lo v are the electron-phonon couplings and 
frequencies, respectively. We have calculated the ener- 
gies of the lowest singlet and triplet states for a free Ceo 
molecule. The results are shown in Fig. |l|. Experimen- 
tally it is found that the low spin state wins for both 
A3C60 and A4C60, A4C60 being a nonmagnetic insula- 
tor 16 and NH3K3C60 being antiferromagnetic with a 
moment (0.7 [Ib per molecule) p£| which corresponds to 
a spin 1/2 system. For A4C60 the triplet-singlet splitting 
is estimated to be 0.1 eV fl6[| . In Fig. [i] this splitting 
is obtained for K ~ 0.07 eV. As expected this value is 
smaller than the value K = 0.11 eV deduced from the HF 
calculation. We observe that a change in the estimate of 
g of course would also lead to a change in this empirical 
number of K. Since, however, the present estimate of K 
seems reasonable relative the the HF value, this calcula- 
tion also gives some support for the value of g deduced 
from PES. 
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FIG. 1. The triplet-singlet splitting of a Cg 4 molecule with 
one H s phonon and an exchange integral K. The arrows show 
the value of K deduced from a HF calculation and the empir- 
ical value deduced from the experimental triplet-singlet split- 
ting. 



IV. A TWO-SITE MODEL 



In this section we study a two-site model to gain un- 
derstanding of a) the different influences of an A g and 
a Jahn- Teller phonon on the metal-insulator transition 
and b) the competition between the Jahn- Teller effect 
and the Hund's rule coupling. We therefore compare the 
case of an A g phonon coupling to a two-fold degenerate 
level (Ax E problem) with the Jahn- Teller E x E prob- 
lem. Thus we study the model 



H 



2 

E 



— phon (i) + H v (i)] + H hop , 



(18) 
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where H c i_ p h on (i) (Eq. (|])) describes the electron- 
phonon interaction on site i, H\j(i) (Eq. (Q)) describes 
the Coulomb interaction on site i, and 

2 

#hop = -t ^ Y^lma^ma + V'LwV'w), (19) 
m— 1 er 

describes the hopping between the two sites. We have 
assumed that there is only hopping between equal m- 
quantum numbers. 

The band gap E g (n) for the filling n is given by 

E g (n) = E(2n - 1) + E{2n + 1) - 2E(2n) (20) 

for a two-site model, where E(N) is the ground-state en- 
ergy of a system with TV electrons. It is useful to consider 
Eg in the limit when U is very large, and then to extrap- 
olate to intermediate values of U. The limit of a large U 
is particular transparent, allowing for simple calculations 
and a qualitative understanding, but it is still relevant for 
the intermediate values of U where the metal-insulator 
transition takes place [Q . The the metal- insulator tran- 
sition happens when E g = 0, i.e., we find the critical 
value of U for which this condition is satisfied. The two- 
site system is much too small for obtaining reliable quan- 
titative estimates of the critical U and the extrapolation 
of the large U results is questionable. Nevertheless, this 
approach can give a qualitative understanding of more 
realistic calculations We consider the limit 

g « lu «W «U 

K^^- = E.j T « W, (21) 
ojo 

where the electron-phonon coupling and the exchange in- 
tegral are just weak perturbations to the hopping and to 
the Coulomb interaction. 

We first study the case of an A g phonon and calculate 
the gap for the filling n = 1. It is then convenient to 
transform the coupling term to the form 

2 

3 ^K-l)(6 l + fe J t ), (22) 
<=i 

where n 4 = J2 mc r Wm^im* and we have neglected an 
irrelevant term — (g 2 /u>o) ^2A2rii — 1). For a system with 
two electrons and in the limit of a large U, hopping is 
almost completely suppressed, and each site has almost 
exactly one electron. The term ( p2| ) then has essentially 
zero coupling to this state, and we obtain 

£ A '(2)«0. (23) 

Since there is only one electron per site, there are no mul- 
tiplet effects. To obtain the energy of a system with one 
electron, it is convenient to transform to bonding and an- 
tibonding operators, e.g., ip± m a = {ipima ± ^2ma)- This 
gives the interaction term 



^[(V'+ma^+mtr + 1pl ma ll>- ma ) (&+ + &+) 

Hi>+m<ri>-m„ + ^ mCT V+ma)(&- + &!)]. (24) 

The bonding level is then occupied by one electron, giving 
the kinetic energy —t. The leading contribution to the 
energy due to the electron-phonon interaction is — \Ejt-, 
giving 

E A °(l) = -t-^E JT . (25) 

In the limit studied here (ujq « W), the phonons react 
to the average electronic charge, which gives rise to the 
factor |. 

We next consider the system with three electrons. The 
wave function for g = and K = is 

|0) = -^d'll^ln^Lt +^11^21X^121 +^211^11^121 
+^211^211^21 + ^211^11^221 + ^IIT^II^t)!™ )- ( 26 ) 

The corresponding energy is (0|£f|0) = — 2t. For a gen- 
eral lattice and for a state which is half-filled apart from 
one electron or one hole, we expect the hopping to be 
enhanced by roughly a factor ^N^cg, where JVd eg is the 
orbital degeneracy JlJ]. However, for the two-site sys- 
tem this enhancement is N^cg fl^| . The effects of the 
electron-phonon coupling are now included in perturba- 
tion theory. For the parameter range considered here 
(u>o << W), there is only coupling to the states fot j 0) . 
Furthermore, the Coulomb energy is given by (0|-Hu|0), 
since there are no states with the kinetic energy -2t which 
couple to |0) via H\j. The energy is therefore 

E A ° (3) = U xx -2t- \e jt - ^-K. (27) 

The corresponding band gap is given in Table |ll| (AxE), 
which also shows the result for the filling n = 2. The ta- 
ble shows that the band gap is reduced by the electron- 
phonon coupling. Extrapolating to intermediate values 
of U then suggests that the electron-phonon coupling in- 
creases the critical value of U where the metal-insulator 
transition takes place. The reason is that the electron- 
phonon interaction in the model (^2|) is not effective for 
integer filling in the large U case, since the charge fluctu- 
ations arc then suppressed. For the states with an extra 
electron or hole, on the other hand, there is a fluctu- 
ating charge due to the hopping of the electron or hole 
and a corresponding lowering of the energy due to the 
electron-phonon coupling. 

The multiplet effects trivially reduce the gap for n = 1, 
since they are effective for the state with three electrons 
but not for the other two states. For n = 2, however, they 
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increase the gap. The reason is that in the integer oc- 
cupation case (N = 4), hopping plays no role for a large 
U and the state can adjust to use the multiplet effects 
optimally. This is not possible for the case of an extra 
hole (N = 3) or an extra electron (N = 5), since then the 
states first of all adjusts to optimize hopping, and only in 
the second place adjust to optimize the multiplet effects. 
The increase of the gap due to the exchange integral has 
been observed earlier |l3| . 

We next consider the Jahn- Teller case of two-fold de- 
generate phonons. We first calculate the energy of the 
state with one electron per site in the large U limit. Per- 
turbation theory shows that each phonon contributes an 
energy Ejt- Since there are two phonons per site and 
two sites, the energy is 

E E <>{2) 



to a reduction of the gap. As K is further increased, the 
system instead tends to go into high spin states, and the 
gap is increased as K is increased. 

It is interesting to compare the results with a different 
approach. We can calculate an effective on-site £/ e //(n) 
for the filling n as 



U eff (n) = E(n + 1) + E{n - 1) - 2E{n), 



(30) 



where E{N) now refers to the energy of a free molecule, 
calculated in Sec. |l[ For K < 2Ejt and for the E x E 
case, we then obtain 



-4E 



JT: 



(28) 



U eff (l) = U av (2)-4E JT + -K 

16 

U eff (2) = U av + 12E JT -—K. 



(31) 



where E ff labels the Jahn- Teller case. Although charge 
fluctuations are completely suppressed in this case, the 
system can still gain energy by introducing a (dynami- 
cal) Jahn- Teller distortion. This is in strong contrast to 
the A g phonons, which only couple to the net charge on 
a given site. In perturbation theory, the energies of the 
states with an extra electron or hole arc 



E E "(1) = -t-E 



JT 



E L 



'(3) = U xx — 2t — Ejt — —K. 



(29) 



The electron-phonon interaction lowers the energy twice 
as much as in the A g case (Eqs. ( |25| ) and (|27|)), simply 
because the phonon is two-fold degenerate and the cou- 
pling therefore is to twice as many phonons as in the A g 



It is then tempting to assume that for the filling n = 
1 and n = 2 the gap is given by U e ff(l) — 3i and 
U e ff(2) — 4t, respectively. Comparison with Table III 
shows that this is incorrect in the limit studied here (Eq. 
(pl|)). In particular, for n = 1, this approach would even 
predict the wrong sign for the the electron- pho non con- 
tribution to the gap. The reason is that Eq. ( |30[) assumes 
that the electron-phonon and Hund's rule couplings can 
adjust to the instantaneous occupation of a given site. 
As we have seen above, this is not possible in the limit 
(Eq. (pl|)) for the states with an extra electron or hole. 
U e ff(2) may instead become relevant in the limit when 
the exchange coupling K and the Jahn- Teller energy Ejt 
are much larger than the hopping energy. 

In a similar way we have calculated the gap for the 
T x H two-site problem. The results are shown in Table 
[Tv| . These results also illustrate the competition between 
the Jahn- Teller effect and the Hund's rule coupling. 



case. The corresponding gap is shown in Table III, which 
also shows the gap for the filling n = 2. 

For n = 2 and K < Ejt the Jahn- Teller effect domi- 
nates over the Hund's rule coupling. In this case as well 
as for n — 1 the gap is then increased by the electron- 
phonon interaction, while in the A g case it is decreased. v - DYNAMICAL MEAN-FIELD CALCULATIONS 
The reason is similar as for the multiplet effects in the A g 
case, discussed in the second paragraph below Eq. (27). 
For the parameter range (Eq. (pl|)) we are considering 
and for integer filling, the hopping is very efficiently sup- 
pressed, and the Jahn- Teller system can therefore adjust 
efficiently to the electron-phonon interaction. With an 
extra electron or hole, however, the system cannot take 
advantage of the electron-phonon interaction to the same 
extent, since the wave function primarily optimizes the 
hopping of the electron or hole. According to Eq. (20), 



this leads to an increase of the gap. In the A g case, on 
the other hand, the system cannot couple to the phonons 
in the integer filling case, and therefore even the reduced 
coupling to the phonons in the case of an extra electron 
or hole is sufficient to reduce the gap. 

Table III furthermore illustrates the competition be- 



tween the Jahn- Teller effect and the Hund's rule cou- 
pling. For K < 2Ejt and n = 2 an increase of K leads 



We now consider the T x H problem in the dynam- 
ical mean- field theory (DMFT). We formulate the infi- 
nite dimensional limit on the Bethe lattice, where the 
nearest-neighbor hopping integral ti m .j m ' is rescaled as 
t* / \fz with the connectivity z going to infinity. We fur- 
ther simplify the hopping by setting U m ,jm' oc S rnm > only 
allowing the diagonal hopping. The unit of energy is cho- 
sen such that the bandwidth, W, is set to 2. The Jahn- 
Teller phonons embedded at each lattice sites of the bath 
can be easily incorporated into the effective medium of 
the impurity Anderson model p9[ , since they are Ein- 
stein phonons without direct intersite couplings between 
them. Therefore the lattice contribution of phonons is 
implicitly included in the medium electron Green's func- 
tion. 

We here focus on the problem where the exchange in- 
tegral K = 0. The effective impurity Anderson model is 
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solved using the quantum Monte Carlo (QMC) technique 
with the Fye-Hirsch algorithm |30|, which has only mild 
"fermion sign-problems" with our Hamiltonian. Here we 
treat fully quantum mechanically the phonon fields which 
are updated together with the fermion auxiliary fields in 
each Monte Carlo step. Details of the implementation of 
the QMC technique in the DMFT can be found elsewhere 
in the literature |l9| , |29| ]. We have used the discretiza- 
tion step for the Trotter breakup, At = 1/3, throughout 
this paper unless mentioned otherwise and more than one 
million Monte Carlo sweeps are taken for each iteration 
of the self-consistency loop. For the case where the ex- 
change integral K > the QMC method has a serious 
"sign-problem", and this case is treated in Sec. [v|D us- 
ing exact diagonalization. 



VI. RESULTS 

We study how the metal-insulator transition depends 
on the parameters, e.g., the ratio U/W for different 
strengths of the electron-phonon coupling A and the ex- 
change integral. When the system becomes insulating a 
gap is opened up in the electron spectral function A(u>). 
This shows up in the electron Green's function G(r) cal- 
culated for imaginary times r. For instance if 



0.12 



A{u) = 



for 



< A, 



we obtain 



G(t = /3/2) < e- A/3 / 2 , 



(32) 



(33) 



where (3 = l/T. Thus G((3/2) decays exponentially with 
(3 for an insulator. We therefore use the behavior of 
G((3/2) as a measure whether the system is a metal or an 
insulator. It is also interesting to study the charge fluc- 
tuation ((n — no) 2 ), which is an average of ((n, — «-o) 2 ) 
and where no is the average occupancy per site. This 
quantity is expected to become small but nonzero at the 
metal-insulator transition. 




FIG. 2. (a) G(/3/2) for p = 16 and (b) the charge fluctua- 
tion ((n — 3) 2 ) as a function of U/W for different values of A 
and coupling to H 9 phonons in AsCeo- 



It is interesting that U c as a function of A initially is 
reduced very strongly as A is increased, while for larger 
values of A the decrease is slower. This can be understood 
in terms of the results in Eq. ( pT|, p"3| ) for the Jahn- Teller 
energy in the weak- and strong-coupling limits. This 
shows that the electron-phonon energy increases much 
faster with A (by a factor of 5/2) in the weak-coupling 
limit than in the strong-coupling limit. This is particu- 
larly relevant for the large U integer filling case, where 
the electron-phonon interaction has a similar effect as in 
the free molecule, and it gives a qualitative explanation 
for the dependence of U c on A. 



A. Jahn- Teller H s phonons 

We first study the Jahn- Teller H 9 phonons and con- 
sider the case of half-filling, i.e., three electrons per site. 
Fig. H shows G(/3/2) as a function of U/W for differ- 
ent values of A. The figure illustrates that for a given A, 
G(/3/2) is reduced as U is increased and at some critical 
value U c , G((3/2) becomes very close to zero (not exactly 
equal to zero due to the finite temperature), where a 
metal insulator transition takes place. The figure also 
illustrates that the charge fluctuations ((n — 3) 2 ) are 
strongly reduced in the insulating state. The critical 
U c /W is reduced as A grows. The reason for this was 
discussed extensively in Sec. [V. 



0.13 



0.12 - 




FIG. 3. G(/3/2) for p = 16 as a function of U/W for 
A = 0.8 and A = and coupling to A 9 phonons in A3C60. 
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B. A g versus Jahn- Teller H 9 phonons 

We next compare the Jahn- Teller H 9 phonons with A g 
phonons. Fig. |^ shows G(/3/2) in a system with cou- 
pling to A g phonons as a function of U for A = and for 
A = 0.8. Comparing the results for A = and A = 0.8, 
we can see that the coupling to A g phonons increases 
the critical U c where the metal-insulator transition takes 
place, while in the case of H g phonons this value is de- 
creased (see Fig. ||). The reason for this change was 
discussed extensively in Sec. [V. 

It is interesting to observe that both the electron- 
phonon interaction and the electron-electron interaction 
by themselves would tend to reduce G(0/2), In a weak- 
coupling theory where we simply add the lowest order 
contribution to the self-energy from each interaction, we 
would then predict that the two interactions work to- 
gether in reducing G(/3/2). This is indeed what is found 
in Fig. U for U/W = 0.5. Not too surprisingly, the argu- 
ments presented below Eq. (|2^) and assuming the large 
U limit are incorrect for U/W = 0.5. It is not surpris- 
ing that these arguments become correct for large values 
of U/W, but it is important that they are qualitatively 
correct already for U/W ~ 1.5, which is on the metallic 
side of the metal-insulator transition. 

For a system with only a Coulomb interaction U, the 
metal-insulator transition can be thought of as resulting 
from the competition between kinetic and potential en- 
ergy. When the coupling to the phonons is introduced 
the hopping of the electrons is reduced, since an electron 
tends to drag a cloud of phonons along its path. It is 
then tempting to assume that the coupling to phonons 
will move the metal-insulator transition towards smaller 
values of U. Our results show that for the parameter 
range considered here, the effect is the opposite for the 
case of coupling to A g phonons. 

The coupling to the A g phonons in A 3 C6o is weak [Q, 
and the A g phonons should not play an important role. 
There is, however, a strong coupling to a charge carrying 
piasmon derived from the t\ u electrons pi] . This pias- 
mon couples in the same way as the A g phonons. Its 
energy (0.5eV) is larger than the phonon energies. How- 
ever, whether we consider the piasmon energy to be small 
or large we arrive at the conclusion that it increases U c . 
For small values of the piasmon energy this follows from 
the study of the A g phonons, and for large energies of 
the piasmon, we can introduce an effective U e ft as in 
Eq. (|30|). The coupling to the piasmon reduces f7 e // 
and therefore the metal-insulator transition happens for 
a larger bare U. 



C. A3C60 versus A4C60 

Fig. I shows G{(3/2) for A 4 C 60 . The figure illustrates 
that U c is smaller at filling four than filling three (see 



Fig. |i|). This is not surprising in view of Eq. (|1 
and Table [V, since the energy is lowered more by the 
electron-phonon coupling for filling four than for filling 
three. According to Eq. ( p0| ) this increases the gap of the 
insulating state more for filling four than for filling three, 
i.e., it makes U c smaller for filling four. As discussed in 
the next section, this result is, however, modified by the 
Jahn- Teller coupling. 



0.2 



0.15 - 



CM 

s 0.1 - 



0.05 




1 

U/W 

FIG. 4. G(/3/2) for /3 = 16 as a function of U/W for 
A4C60 and different values of A. 



D. Competition between Jahn- Teller and Hund's 
rule coupling 

We finally discuss the competition between the Jahn- 
Teller effect and the Hund's rule coupling. Since there is 
a sign-problem in the dynamical mean-filed calculations 
if we use the full multiplet coupling in Eq. (Q), we use 
an exact diagonalization technique. This requires that 
the system size is small. Therefore we consider a system 
with just four sites and we consider the E x E problem 
(Eq. (||)). We furthermore limit the Hilbert space by not 
allowing more than two phonons per site. To reduce the 
discreteness of the one-particle spectrum for such a small 
system, we pick the hopping integrals randomly. We cal- 
culate the gap according to Eq. (^0|) and add some finite 
size corrections [Ol 



E r „ e (U xx ) — Eg 



Uav(2) 



E. 



(34) 



where U av {2) / N S i te is a contribution to the gap from the 
electrostatic energy of a finite system, N s ite is the num- 
ber of sites and Eg 7 ^ is the gap of a system without 
any Coulomb interaction. Both these corrections go to 
zero for a large system. Assuming that that E r g ed grows 
linearly with U xx , we can then use 



Uxx E„ [Uxx] 



(35) 



as a crude estimate of the critical U xx . This quantity is 
shown in Fig. |^. The results can be qualitatively under- 
stood from the results in Table III for a two-site system, 
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although the parameter range considered here is outside 
the range where the results in the table are valid. For 
A = 0, the critical value of U xx is reduced as K is in- 
creased, as expected However, as A is increased, the 
competition between the Jahn- Teller and Hund's rule ef- 
fects leads to an increase of the critical U xx . This is in 



agreement with Table [II, although the increase is faster 
than for the parameter range of this table. For a par- 
ticular value of A, the critical U xx becomes comparable 
to the value for K = A = 0. As A is further increased, 
the Jahn- Teller effect becomes dominating and the crit- 
ical value of U X x decreases again. This critical value is, 
however, larger than for K — 0, due to the competition 
between the Jahn- Teller and Hund's rule coupling. 




0.25 



FIG. 5. The simple estimate U xx — E r g ed of the critical U xx 
as a function of A for different values of the exchange integral 
K. The one-particle band width is IT = 2 and the phonon 
energy is ljq = 0.5. For A we use A = 4</ 2 / (u)qW) . 



VII. DISCUSSION 

It is interesting to discuss these results in the context 
of A„C6o- Considering just the Hubbard U (g = K = 0), 
it is found that the metal- insulator transition in A 3 C6o 
takes place at the upper range of what is believed to be 
physical values of U/W, while for A4C60 this happens at 
the lower range of these parameters. This agrees nicely 
with the fact that A3 Ceo is a metal but A4C60 is an 
insulator. However, to explain why A4C60 is not antifer- 
romagnetic, we need to include the coupling to the Jahn- 
Teller H g phonons. This substantially lowers the critical 
U c for both systems, and it puts U c /W of A3C60 at the 
lower range of the physical parameters. This puts our 
understanding of A3C60 being a metal into question. We 
find, however, that the competition between the Jahn- 
Teller effect and the Hund's rule coupling increases U c 
again. The coupling to the t\ u plasmons in A3C60 should 
lead to an additional increase of U c in this system. This 



makes it understandable that A3C60 can be a metal 
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TABLE I. The ground-state energy E(N) for N electrons 
in the E x E case. The quantity N(N - l)J7 „(2)/2 has been 
subtracted, where U a v{2) = U xy + K/3 is the average inter- 
action for a full shell. 



N 


E(N) = E(N) 


-N(N - l)£/„„(2)/2 




Low spin 


High spin 


1 


—2Ejt 




2 


-8Ej T + Ik 


— ±K 

3 n 



3 —2E,jT 



TABLE II. The ground-state energy E(N) for N electrons 
in the T x H case. The quantity N(N - l)J7 „(3)/2 has been 
subtracted, where U av {3) = U xy is the average interaction for 
a full shell. 



N 


E(N) = E(N) 


-N(N-l)U av (3)/2 




Low spin 


High spin 


1 


— t;Ejt 




2 


-WEj T + 4:K 


— \Ejt — K 


3 


-^E JT + 2K 


-3K 


4 


-WE JT +4K 


— \Ejt — K 


5 


-\E JT 





TABLE III. E g (n) - U av (2) - d 2 (n)t for the A x E and the 
E x E two-site models as a function of the filling n, where E g 
is the band gap, U av (2) is the average Coulomb interaction 
and (fc(n)t is the hopping contribution, with (fe(n) = -3 and 
-4 for n — 1 and 2, respectively. The results are symmetric 
around half-filling n = 2. 



Syst 


n 


E g {n) -U av {2) 


- d 2 (n) 






K < 2Ejt 


K > 2Ejt 


Ax E 


1 


— Ejt 




Ax E 


2 


—Ejt + ^ 


K 


E x E 


1 


QEjt 




ExE 


2 


30Ej T -™K 


2Ejt + fK 



TABLE IV. E g {n) - U a v{3) - d 3 {n)t for the T x H two-site 
model as a function of the filling n. The hopping contribu- 
tion to the gap is given by d(n)t, where ds(n) — -3, -5 and -6 
for n = 1, 2 and 3, respectively. The results are symmetric 
around half-filling n — 3. 



E g (n) -Uav{3) -d 3 {n)t 


n 




Ejt %Ejt < K < \E JT 


K > \E JT 


1 




5E JT + 4K 


fEjT - K 


2 
3 


35S./T - 
55 rn 

— Ejt 


- fK 5E JT + fK fEjT + 3K 
-8K -\e jt + 12K 
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